Single-Shot Readout of Macroscopic Quantum Superposition State in a 

Superconducting Flux Qubit 
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Single-shot readout experiments were performed on the two lowest-energy states of a supercon- 
ducting qubit with three Josephson junctions embedded in a superconducting loop. We measured the 
qubit state via switching current (7sw) of a current-biased dc-SQUID, a quantum detector surround- 
ing the qubit loop. The qubit signals were measured in a small Jsw regime of the SQUID, typically 
less than 100 nA, where the Jsw distribution is particularly narrow. The obtained single-shot data 
indicate that the qubit state is readout, through the flux generated by the qubit persistent-current, 
as energy eigenstates rather than current eigenstates. 

PACS numbers: 03.67.Lx, 85.25.Dq, 85.25.Cp, 03.65.Yz 



Among many candidates toward realization of quan- 
tum computation, superconducting qubits based on 
Josephson Junctions take on a position of increasing im- 
portance Although their coherence time still need 
to be substantially improved, solid state qubits such 
as semiconductor qubits or superconductor qubits have 
the distinct merit of scalability due to mature state of 
nanometer-scale fabrication technology. Recently, sev- 
eral groups have achieved a fairly long coherence time, 
and gate operation for single qubits has become available 
by the application of custom made resonant microwave 
pulses, which control Rabi oscillations 2, 40ilEIE0- 
Rotational gate operation together with two-qubit con- 
trolled NOT operation has also been successfully demon- 
strated 111- As regards the qubit readout, averaging pro- 
cesses are commonly used and a way to overcome imper- 
fect readout visibility has been a topic of debate. 

In this paper, we report the observation of the single- 
shot readout of a superconducting flux qubit com- 
prising three Josephson junctions in a loop 0. Al- 
though, single qubit operation and entanglement between 
two qubits have been demonstrated in superconducting 
charge qubits, fabricated Josephson circuits exhibit static 
and dynamic charge noise due to background charged 
impurities. In contrast, in flux qubits there is relatively 
little magnetic background noise. Superconducting flux 
qubits therefore have the potential advantages of a longer 
coherence time and greater stability. 

Figure 1(a) shows a scanning electron microscope im- 
age of a sample; a superconducting qubit (inside loop 
with three Josephson junctions of critical current /q"t>it ^ 
0.6/xA for larger junctions) and an under-damped dc- 
SQUID (outside loop with two small Josephson junctions 
of critical current /f*^ ~ 0.15/iA for each junction) as a 
detector. The qubit contains three Josephson junctions, 
two of which have the same Josephson coupling energy 
Ej = hlcl2e. The third has aEj, with 0.5 < a < 1. 
The a value can be controlled by the ratio of the area of 
the smallest Josephson junction to the other two larger 
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FIG. 1: (a) Scanning electron microscope image of a supercon- 
ducting qubit (inside loop) and a dc-SQUID (outside loop) as a 
detector. Arrows indicate the directions of the supercurrent in the 
qubit. (b) Bird's-eye view of the Josephson energy of the qubit in 
phase space, (c) Schematic drawings of the measurement circuit. 
The squares with crosses represent Josephson junctions. 



junctions in the qubit. We used a sample with a — 0.8 
with the areas of the larger and smaller junctions given 
by 0.1x0.4 fim^ and 0.1x0.32 ^m^, respectively. The 
loop area ratio of the qubit to the SQUID was about 
3:5. Arrows represent macroscopic supercurrent flows in 
the qubit. The dc-SQUID with two Josephson junctions 
was under-damped (no shunt resistor). The qubit and 
the dc-SQUID were coupled magnetically via mutual in- 
ductance M ~ 7 pH. The aluminum Josephson junctions 
were fabricated by suspended bridge and shadow evapo- 
ration techniques liOj. The sample was cooled to 25 mK 
with a dilution refrigerator and it underwent a supercon- 
ducting transition at ~1.2 K. In order to reduce external 
magnetic field fluctuations, both the sample holder and 
the operating magnet were mounted inside a three-fold /.t- 
metal can in the refrigerator. Figure 1(b) shows the qubit 
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FIG. 2: (a) Switching current of dc-SQUID detector. / = 
*I'oxt/4"o, where "I>oxt is the external magnetic flux penetrating the 
qubit loop.(b) Schematic illustration of the qubit energy dispersion 
relation, (c) Expected readout for the qubit energy eigenstates |0) 
and (d) Enlarged view near / = 0.5. (e) Enlarged view near 
/ = 1.5. The two states of the qubit are read out distinctly. 



FIG. 3: DC-SQUID switching current (/sw) distribution as a 
function of /sq = 'I'sq/'I'o, where 'I'sq is the external magnetic 
flux penetrating the SQUID loop. The standard deviation of the 
measured switching-current IsD (•) a^nd its ratio to /gw; Isd/Ibv^ 
(o) are plotted. Note that I^d suddenly drops in regions where 
/s„ is smaller than ~100 nA. 



potential energy [/(7-|_,7_) in the phase space, where 
-— ~ 2 + a — cos 71 — cos 72 — acos(27r/ — 71 — 72), (1) 

and 71+72 + 73 = 27r/, with 7,- being the phase difference 
across the three junctions [ij . Junctions 1 and 2 have 
identical iJj, but junction 3 with 73 has aEj. The filling 
/ is the external flux $ext penetrating the qubit loop 
normalized by the superconducting flux quantum $0 = 
In a low-energy approximation, the qubit potential 
can be described by 7+ alone, where 7± = (71 ± 72) /2. 

By carefully designing the junction parameters [211131 j 
the inner loop can be made to behave as an effective two- 
state system In fact, the read-out result of the 
qubit changes greatly with qubit design, r ang ing from the 
purely classical to the quantum regime The qubit 
is described by the Hamiltonian Hq = ^{sfaz + ^ctx), 
where ax.z are the Pauli matrices. Two eigenstates of (t^ 
are localized states with the supercurrent circulating in 
opposite directions, i.e., the clockwise state \R) and the 
counter-clockwise state \L). The dc-SQUID picks up a 
signal that is proportional to tr^ . The barrier for quantum 
tunneling between these states depends strongly on the 
Josephson energy Ej, charging energy Ec and a values. 
The tunneling matrix element is given as 

tLR^^^ 1.3/^ exp(-0.64W ^) (2) 
z V 

for a = 0.8. In the present sample, Ej — 1.27 meV 
and E_j/Ef. — 50. As schematically shown in Fig. 1(c), 
the qubit is biased with a static magnetic flux $ext us- 
ing an external coil. A microwave antenna was placed 



just above the sample chip to induce oscillating mag- 
netic fields in the qubit loop. The switching voltage 
of the SQUID was measured by the four-probe method. 
The detector SQUID was current biased through resis- 
tors Ri to avoid the parasitic capacitance of the leads. 
The measurement lines were composed of a pair of re- 
sistors (i?/ — Rv ~ 200ri), constantan twisted wires up 
to the mixing chamber position, fiexible coaxial cables 
from the mixing chamber up to the top flange and 10 nF 
through capacitors at room temperature. The SQUID 
bias current (/b) in a saw-tooth wave was applied with a 
repetition rate of few hundred Hz. The amplitude of /b, 
typically up to a few hundred nA, depends on the qubit 
operating point. The SQUID switching voltage was am- 
plifled by a differential amplifier and (/gw, Vsw) switching 
was recorded when a voltage greater than the threshold 
{'-^ 30/iV) was first detected during every Jb ramp. 

Figure 2(a) shows the modulation of the measured dc- 
SQUID switching current against the external magnetic 
flux as a function of / = <I>ext / 'I'o ■ Each dot in the graph 
corresponds to a single-shot readout without averaging. 
The classically stable current eigenstates \L) and \R) are 
no longer stable in the presence of quantum tunneling 
Ilr. As shown in Fig. 2(b), the qubit ground state |0) 
and the first excited state |1) change the direction of the 
persistent current as a function of /. Here, 

|0) = sin(^)|L)+cos(^)|ii') (3) 

|1) = -cos(^)|L)+sin(^)|i?), (4) 

where Of — arctan(A/ey). Thus, the energy eigenstates 
|0) and |1) are superpositions of macroscopically distinct 
states \L) and \R) 0,0. The expected qubit response 
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FIG. 4: (a) Normalized qubit steps as a function of / = 'I'cxt/'3?0) 
where <I>ext is the external magnetic flux penetrating the qubit loop. 
Temperature-averaged expected qubit responses are also shown as 
solid curves, from which the sinusoidal background modulation of 
the SQUID has been subtracted. The step height is normalized 
at 2. Each curve is shifted vertically for clarity, (b) The classical 
behavior of a two-level system observed in a different sample. 



in Oz is shown schematically in Fig. 2(c). Therefore, 
this step-like feature, which is called a "qubit step", is 
expected to occur near / « i -I- n, where n is an inte- 
ger (indicated by the boxes in Fig. 2(a)). Figure 2(d) 
shows an enlarged view near / = 0.5, where the switch- 
ing current distribution of the SQUID is very broad. It is 
almost impossible to distinguish between the two states 
of the qubit |0} and |1) by a single-shot measurement. 
In Fig. 2(e) near / = 1.5, the switching current distri- 
bution of the SQUID is narrow enough for us to observe 
the x-shaped qubit steps clearly. The observed amount 
of magnetic-flux shift A$ shown in Fig. 2(e) is consis- 
tent with the expected value 2aM/;?"''" = 3.4 x lO^^^o- 
Here, we are able to observe the two energy eigenstates 
of the qubit with a single-shot measurement. 

Figure 3 shows the dc-SQUID switching current dis- 
tribution together with the average switching current as 
a function of /sq — <i>sQ/*i'oj where $sQ is the external 
magnetic flux penetrating the SQUID loop. The sam- 
ple used here is nominally identical to that of Fig. 2 
and shows a similar ^-shaped crossing. We adopted the 
standard deviation as a measure of the switching current 
distribution. We obtained a typical external magnetic 
flux dependence /gw cx | cos ^^^^j. For larger /gw values, 
the standard deviation of the measured switching-current 
isD (•) remains large 4 nA) due to stochastic switch- 
ing processes, which can be described by the macroscopic 
quantum tunneling (MQT) of the "phase-ball" in the 
tilted Josephson washboard potential. The flux detection 



FIG. 5: Normalized dc-SQUID switching current readout (open 
circles) after taking an average at each flux bias. The solid curves 
represent thermal-averaged theoretical values. 

performance of the SQUID is determined by the distribu- 
tion and sensitivity of the detector dc-SQUID. The flux 
sensitivity of the SQUID is proportional to the slope of 
the /sw modulation curve, i.e., | J^^l and becomes higher 

near /sq « ^ -f to, where m is an integer. In addition, 
the switching current distribution becomes particularly 
narrow in the small critical current region, where /sw is 
typically smaller than 100 nA. This region corresponds 
exactly to the crossover region E]/Ec ~ 1 of the detec- 
tor dc-SQUID, since the area of each Josephson junc- 
tion in the dc-SQUID is 0.1/^mxO.l/im, which is roughly 
four times smaller than that of the qubit junction, re- 
sulting in Et,IEc ~ 50/4^ ~ 3 near /sq ~ 0, where 
/sw — 300 nA; therefore Ej/Ec < 1 when Igw is smaller 
than ^ 100 nA. It is notable that a precipitous change 
in the signal-to-noise ratio occurs at Ej/Ec ~ 1. An- 
other possibility is dissipation induced narrowing. The 
ground state wavefunction of the SQUID localizes more 
strongly by the incoming noise, resulting in a narrower 
/sw distribution^^!- For small-capacitance low-critical 
current Josephson junctions, it is known that the width 
of switching histogram decreses as a consequence of acti- 
vation process due to effective temperature over the dissi- 
pation barrier [l9l |. Consequently, we have chosen to use 
this extremely sensitive region of the SQUID detector. 
By choosing appropriate integers m and n when design- 
ing the qubit and SQUID, we succeeded in tuning the 
qubit operating point, / — 1.5, to the high sensitivity 
region of the SQUID. 

By subtracting the background modulation of the dc- 
SQUID, we can extract the signals from the qubit as 
shown in Fig. 4(a). The normalized switching current 
of the dc-SQUID against the external magnetic flux is 
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shown as a function of / = $cxt / at different temper- 
atures. Each dot in the figure corresponds to a single 
measurement and no averaging is performed. The sohd 
curves are fitted to the ground and excited states with 
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where e/ = 2/p$o(/ - 1.5) with /p ~ 4''"'^'^ 1 - 
The thicker stripe corresponds to the ground state |0) 
and the other stripe corresponds to the excited state 
Here, we should point out that we can expect the dc- 
SQUID to detect \L) or \R) , i.e., the current eigenstate, 
by a single-shot readout, in the same way as the Stern- 
Gerlach apparatus. However, the state observed by this 
single-shot measurement is the energy eigenstate not the 
current eigenstate of the qubit. Recently, we have the- 
oretically investigated the switching current behavior as 
a function of the external magnetic flux <i>aubit/'I'o with 
the conditions of our measurement method 14]. We found 
that the qubit energy eigenstate of superposition between 
|L) and |R) is maintained until the readout, i.e., SQUID 
switching event, if the relaxation time is long enough 
and qubit-SQUID interaction is weak as is in our experi- 
ments. Whereas, in the superconducting charge qubit, 
single-shot detection on a non-energy eigenstate basis 
has recently been reported |22|. We recorded two thou- 
sand switching events under each fixed external magnetic 
field. The measurement results emerge randomly on ei- 
ther stripe and there is only a small fraction of points 
in between. Several events were detected in the excited 
states at the lowest temperatures, which may be due to 
unexpected excitation related to certain cavity modes. 

For comparison. Fig. 4(b) shows an example of a classi- 
cal two-level system observed in a different sample. The 
sample parameters are Ej/Ec = 185, Ej — 2.4 meV, 
a — 0.8, with a calculated qubit energy gap A = 0.2 ^eV. 
In this sample, the readouts for the two states are com- 
pletely separated as \L) or \R) by a single shot measure- 
ment. Moreover, the readout does not contain x shaped 
crossing qubit steps. This means that the tunneling term 
A of the qubit is so small that decoherence easily destroys 
the coherent tunneling between the wells. As the MQT 
is suppressed by the high tunnel barrier, only thermal 
excitation or incoherent tunneling is possible and, conse- 
quently, classical bi-stable branches are observed. 

If we take the average at each flux bias /, the steep 
step observed at lower temperatures becomes rounded at 
higher temperatures as expected. Figure 5 shows, the 
normalized readout after averaging (open circles). This 
averaging corresponds to time ensemble averaging. Ther- 
mally averaged expectation values shown as solid curves 
in Fig. 5 are calculated from the canonical distribution : 






■>I1 


RF 3 B5 GHi 


















|0 













50 100 150 200 250 

count 






L..i.i.>. 








50. 

58' 







RF off 




57' 






|0 




56' 











50 100 1 


2[ 


25 



FIG. 6: A precise readout distribution measurement of the region 
indicated in Fig. 4(a) (25 mK) (a) with and (b) without microwave 
irradiation. This sample is from the same set of samples showing 
the x-shaped crossing used in Fig. 4(a). The occupancy of the 
excited state, |1), increases with a microwave at a resonant point. 
The corresponding histograms at the resonance point, / = 1.50102, 
are shown in the right panels. 
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where £'o(i) = -(+)i^e/2 + A2, and = ksT is the 
effective temperature, fee is the Boltzmann constant. We 
obtained reasonable agreement between the experimen- 
tal and calculated results. This is an example of the 
direct observation of the thermal distribution of a coher- 
ent macroscopic quantum object. This shows that even 
a macroscopic quantum object obeys canonical distribu- 
tion as if it were a microscopic object. We observed two 
small dips and two small peaks in the averaged qubit step 
especially at lower temperatures. These dips and peaks 
appeared in a symmetrical position and could be due to 
unexpected resonant modes in the environment formed 
by the sample and sample holder. 

Figure 6 shows a density plot and a histogram of the 
switching current of the sample showing the x cross with 
and without a microwave. The left panel shows an en- 
larged view of the region indicated by the box in Fig. 
4(a) at 25 mK, and the right panel shows a histogram at 
the resonance point, / — 1.50102. The readout is clearly 
separated into |0) and |1) states even with a continuous 
microwave. Without a microwave the readout is largely 
in the |0) state. The residual probability in the |1) state 
is due to incoherent thermal excitations and is very small. 
In contrast, we found that the probability appeared in the 
|1) state with microwave irradiation. Ideally, the proba- 
bilities of detecting these two states should be the same. 
A microwave transition enables us to estimate the coher- 
ence time in the system. Tg* — 5 ns has been obtained 
from the line width of the transition under a resonant 
microwave. The sample showing no x cross (Fig. 4(b)) 
did not respond to the microwave at all, while a coherent 
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transition by the resonant microwave was observed for 
the sample showing x cross (Fig. 4(a)). The observa- 
tion of the transition induced by the microwave indicates 
the existence of macroscopic quantum coherence in the 
system. This is a consequence of a tunnehng term in 
the Hamiltonian. Therefore, in addition to the x cross- 
ing qubit step, the observation of a microwave resonant 
transition provides additional evidence for the coherent 
superposition of macroscopic supercurrent states \L) and 
\R). 

In conclusion, we reported on the single-shot measure- 
ment of a macroscopic quantum state in a superconduc- 
tor qubit. We achieved a high-resolution readout by care- 
fully selecting the flux bias point and taking advantage of 



the highly sensitive switching current region of the detec- 
tor dc-SQUID. We believe that single quantum event de- 
tection in fabricated quantum circuits will contribute to 
a further understanding of the quantum/classical bound- 
ary including quantum measurement processes, powered 
by the recent rapid development of nano-technology. 
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